Spectral and Transport Properties of D-Wave Superconductors With Strong 

Impurities 
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One of the remarkable features of disordered d-wave superconductors is strong sensitivity of 
long range properties to the microscopic realization of the disorder potential. Particularly rich 
phenomenology is observed for the - experimentally relevant - case of dilute distributions of isolated 
impurity centers. Building on earlier diagrammatic analyses, the present paper derives and analyses 
a low energy effective field theory of this system. Specifically, the results of previous diagrammatic 
T-matrix approaches are extended into the perturbatively inaccessible low energy regimes, and the 
long range (thermal) transport behaviour of the system is discussed. It turns out that in the extreme 
case of a half-filled tight binding band and infinitely strong impurities (impurities at the unitary 
limit), the system is in a delocalized phase. 
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I. INTRODUCTION 



Disorder in the cuprates is commonly attributed to the 
presence of dilute distributions of strong local scatter- 
ers - e.g. Zn-atoms replacing Cu - immersed into the 
two-dimensional matrices of Cu02 planes Q|. While for 
conventional superconductors details of the microscopic 
realization of the disorder potential are of little, if any 
concern, the situation with anomalous superconductors 
is different. In fact, microscopic aspects of the impurity 
scattering have proven responsible for many of the noto- 
rious difficulties encountered in developing quantitative 
theories of disordered quasi-particles in d-wave supercon- 
ductors: 

First, the spectral characteristics of d-wave super- 
conductors polluted by dilute strong scatterers - com- 
monly dubbed 'unitary scatterers' - differ strikingly from 
those of systems with continuously distributed disor- 
der: While in the latter case, the spectral density of 
quasi-particles categorically vanishes upon approaching 
zero energy ^,||, it may diverge for unitary disorder |Q. 
This is demonstration of the fact d|| that the stan- 
dard paradigm of "insensitivity of long range properties 
of fermion systems to details of the disorder" is essen- 
tially violated in anomalous superconductors. 

Second, conventional diagrammatic perturbation the- 
ory rules out as a tool to quantitatively address the low- 
energy regime of the system. Based on self-consistent 
T-matrix (SCTA) approximation schemes, early pertur- 
bative approaches to the problem |t]|| did not take the 
existence of infrared-singular diagram classes J9|,[l0| into 
account. These diagrams, similar in appearance to the 
diffusons and cooperons of normal metals, reflect the ex- 
istence of a Goldstone mode linked to the particle-hole 
exchange symmetry of the system. For not too small en- 
ergies, they can be brought under control by summing up 
weak localization type corrections to the Gorkov Green 
functions ||,|l0|]. However, for low energies these correc- 
tions proliferate and perturbation theory produces un- 



physically divergent results. For the same reason, dia- 
grammatic theory cannot be applied to assess the trans- 
port behaviour of the quasi-particle system. 

In contrast, field theoretical approaches to the prob- 
lem ^j,||j5| provide the means for controlled inclusion of 
all low-energy soft modes. However, third, the specific 
case of dilute impurity distributions, commonly referred 
to as Poisson distributed disorder, has been out of reach 
of previous fieldtheoretical formulations. (This is because 
the standard construction route of field theories of disor- 
dered fermions is based on models of continuously, mostly 
Gaussian distributed disorder.) 

To surmount these difficulties, a number of alter- 
native approaches to the problem have been devel- 
oped. Building on an explicit representation of the 
realization-specific Gorkov Green function, Pepin and 
Lee (PL) |Q found the low energy single particle den- 
sity of states (DoS) v(e) to be power-law divergent, 
v{e) ~ l/|eln 2 (|e|)|. While this result was obtained for 
the specific case of a half filled band, PL argue that the 
band center singularity should be insensitive to not too 
strong variations of the chemical potential. (Considering 
applications to the cuprates, the stability issue is of prime 
relevance: At half filling, or zero doping, the Cuprates are 
in an anti-ferromagnetic phase, i.e. the very formation 
of a superconducting phase requires a finite amount of 
dopands. This implies that strong sensitivity of the band 
center singularity to deviations off half filling would im- 
pose severe limitations on the practical relevance of the 
result.) 

However, this latter statement is at variance with nu- 
merical Q and analytical [^0) findings indicating that 
the band center singularity represents a highly fragile 
phenomenon. As shown by Yashenkin et al. (YAGHK) 
|I0| the appearance of a zero energy divergence is tightly 
linked to the formation of a second class of soft modes, 
generated by the nesting symmetry of the half filled 
band and existing in parallel to the standard particle- 
hole modes. Within perturbation theory, these modes 
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lead to the formation of positive logarithmic (but not 
power law, like in PL) corrections to the low energy 
DoS. Further, even slight deviations from the combined 
limit (unitarity/half filling) generate a mass gap for the 
anomalous soft modes and induce a crossover to the neg- 
ative logarithmic corrections characteristic for the stan- 
dard particle- hole modes. 

While the perturbative analysis reveals the existence 
of an anomalous phase in the half filled system, it can- 
not reliably predict the low-energy profile of the DoS, the 
critical properties of the anomalous phase, and its trans- 
port properties. To shed more light on the low energy 
phase of the Poisson disordered system in general, Cha- 
mon and Mudry suggested a duality transformation 
relating the limit of unitary scattering to the opposite 
extreme of Gaussian, or Born disorder. The structure of 
this transformation was suggested on the basis of vari- 
ational equations controlling the high-energy regime of 
the system. However, since the notorious soft modes do 
not at all enter the analysis of Chamon and Mudry, the 
extensibility of the transformation into the fluctuation 
dominated low energy regime remains speculative. 

Building on these previous analyses, it is the objective 
of the present paper to describe the Poisson-disordered 
(i-wave superconductor by fieldtheoretical methods. The 
generalized formalism will then be applied to explore 
spectral and localization properties of the system, at and 
away from half filling. We begin our analysis with a 
brief review of the general field theory approach to the 
e?-wave superconductor, a qualitative discussion of the 
extension to pointlike disorder, and a summary of the 
results the generalized model predicts for spectral and 
transport properties of the system (section |llj) . The next 
two sections contain technical details of the construction, 
first for generic pointlike disorder (section |Hl| ) , then for 
the unitary/half-filled limit (section IV). We conclude in 
section [v| 



II. QUALITATIVE DISCUSSION AND RESULTS 

As with normal systems, the properties of disordered 
superconductors are to a large extend determined by 
symmetries. Specifically, a time-reversal invariant spin- 
singlet superconductor obeys the two relations, 



C:H = -a 2 H T a 2 , 
T : H — H T , 



(1) 



where H is the Gorkov Hamilton operator and <7j are 
Pauli-matrices acting in particle-hole space. The first of 
these relations expresses the general particle-hole conju- 
gation symmetry of the spin-singlet superconductor while 
the second is time reversal. In some sense the symmetry 
class identified through ([j]), termed class CI in Ref. ||, 
plays a role analogous to the standard orthogonal sym- 
metry class of disordered normal systems. 



For a conventional superconductor, the presence of a 
large quasi-particle excitation gap implies that the sym- 
metries (0) only marginally affect low-energy, or long- 
range system properties. However, in a gapless environ- 
ment, as realized, e.g., in the ci-wave superconductor, the 
situation is different. Here, the symmetry C entails the 
formation of long-lived particle- hole excitations: A par- 
ticle and a hole with relative energy difference e interfere 
to form a slowly fluctuating composite mode n(e). As 
e — > 0, the coherence between the particle and the hole 
amplitude becomes perfect, implying infrared singular, 
or massless behaviour of n(e). 

The manifestation of these modes in observable low- 
energy properties can efficiently be explored by field in- 
tegral methods. A fermion-replica field theory describing 
disordered superconductors of class CI was first derived 
by Oppermann [|l2[|l3j l and later rediscovered, then within 
the context of <i-wave superconductivity, by Senthil and 
Fisher || . Supersymmetric implementations of the the- 
ory, needed to explore the perturbatively inaccessible in- 
frared region, have been introduced in Refs. |^|. 

All these models were derived for Gaussian distributed 
disorder. By manipulations similar to those employed in 
the construction of the standard cr-models for systems 
of Wigner-Dyson symmetry, the microscopic model was 
mapped onto an action functional 



Z = / VTexp{-S[T}), 



(2) 



where T is a field taking values in the group Sp(2r) (r: 
number of replicas) for the replica-implementations, and 
in the supergroup OSp(2|2) for supersymmetry. (For a 
precise definition of this group, see Eqs. ( |2^ ) and ( |27| ) 
below.) 

As with the standard cr-models, the structure of the ac- 
tion of this field theory is fixed by a few macroscopic sys- 
tem parameters. Specifically, for a 2d gapless supercon- 
ductor with diffusive quasi-particle dynamics, the (bare) 
action has the form 

S[T] = (3) 

= IT / d2r str ( dTdT ~^ _ 2ie str ( T + T_1 )] ' 

where D is the diffusion constant and vq the DoS com- 
puted in the self-consistent Born approximation (SCBA). 
Importantly, the symmetry of the fields T - which in turn 
is fixed by the physical symmetries (Q) - does not permit 
other contributions to the low energy action. In other 
words, the action is completely fixed by the value of the 
two constants D and v§. Notice that for the d-wave su- 
perconductor, and on the level of the bare theory, 
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where A and t are the order parameter and hopping 
strength, respectively and g s is the spin conductance. I.e. 
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the bare, or Drude spin conductance does not depend on 
disorder concentration |J. 

Starting from the action (^), spectral and transport 
properties of the quasi-particlcs can conveniently be ex- 
plored. However, before turning to the phenomenology 
of the system let us briefly discuss what changes with 
the theory (Q) if the disorder is Poissonian distributed. 
Certainly, the fundamental symmetries ([!]) are insensi- 
tive to the specifics of the disorder distribution. Further, 
the quasi-particle dynamics on intermediate time scales, 
larger than the scattering time but smaller than any scale 
related to quasi-particle localization, will continue to be 
diffusive. This leaves us with two principal possibilities: 
Either there are some fundamental reasons that exclude 
the existence of a local effective action of the Poisson 
disordered system. Or, a low-energy action exists and is 
bound to have the same structure as the action (||). 

The result of a technical derivation detailed in section 
[II will be that the second scenario is realized, i.e. that 
the low-energy phase of the Poisson disordered system is 
again described by . The only difference to the Gaus- 
sian case is that the value of the constant vq is now set 
by the SCTA rather than by the SC-BA density of states. 

That both Gaussian and the Poissonian disorder map 
onto the field theory (|J) implies a number of predictions 
on the long range behaviour of the system: 



> First, the quasi-particle DoS obtains as 
J ,(e) = |Re(str((T + T- 1 )P)) T 



(5) 



where (...) = f VTe~ s ^ stands for functional av- 
eraging and V is a constant projector matrix p4[ . 
In the most basic approximation to the functional 
integral, fluctuations of the fields T are neglected 
and one obtains v{e) = vq in agreement with the 
SCBA/SCTA diagrammatic result. Going beyond 
this approximation, 

> a one-loop integration over fluctuations of the fields 
T leads to a negative correction term || , 



5u(e) 
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In A, 



(0) 



where A 



is an ultraviolett cutoff and a the 



lattice spacing. Eq. flm is in agreement with the 
results obtained by extended diagrammatic series 
summations |fTof , i.e. analyses accounting for the 
presence of singular ladder diagrams. Notice that 
the magnitude of the quantum correction depends 
on the (disorder independent!) scale g s (t/A): The 
smaller the anisotropy t/A > 1, the more impor- 
tant become quantum fluctuations. As concerns 
the 'mass parameter' v^t truncating the logarithm, 
there is an important difference between the Gaus- 
sian and the Poissonian case. While for Gaus- 
sian disorder, vq increases quadratically with the 



strength of the impurity potential, v, it decreases 
as v -1 for strong pointlike scatterers. (For quanti- 
tative expressions, see section III below.), i.e. the 
small parameters of the perturbative approaches to 
the two limits scale reciprocally with the impurity 
strength, which is in agreement with the 'duality' 
hypothesis of Ref . . 

> For systems of finite size L and energies smaller 
than the Thouless energy e = spatial fluc- 
tuations of the fields T are effectively frozen out. 
As with normal systems, a non-perturbative inte- 
gration over the spatially constant zero mode To 
obtains results compatible with RMT behaviour. 
Specifically, the low-energy spectral density scales 
as 



Ke) 



if! 
A' 



(7) 



where A = (vqL 2 ) 1 is the single particle level 
spacing. 

t> Turning to transport, it has been shown [|| that 
under renormalization, the conductance scales from 
its disorder-indepedent bare value to zero, with a (3- 
function (3 — l/(47r 3 g 2 ). This implies that for large 
systems quasi-particles should be localized on the 
scale of some localization length £ ~ exp(47r 2 <7 s ); 
In the thermodynamic limit the system is a (spin) 
insulator. Similarly, 

> in the limit L 3> £, mechanisms of quantum lo- 
calization and spectral repulsion conspire to man- 
ufacture a spectral gap at zero energy: Consider a 
quasi-particle at energies e < D/£ 2 , corresponding 
to time scales larger than the time needed to diffu- 
sively explore a localization volume ~ £ 2 . Confined 
to stay within the localization volume, the quasi- 
particle will exhibit ergodic behaviour, i.e. the dy- 
namics will become random matrix like. One thus 
expects H that 



i/(e) 



J! 

A f 



e < 



D 



(8) 



where Aj = (z4)£ 2 ) -1 is the level spacing of a 
single localization volume. (While for the two- 
dimensional superconductor Eq. (H) has the status 
of a conjecture, similar predictions for quasi one- 
dimensional quantum wires of non-Wigner-Dyson 
symmetry were confirmed by transfer matrix meth- 
ods @@.) 

Summarizing, we find that for both Gaussian and Pois- 
son disorder the system scales to a spin-insulator phase. 
While the DoS's obtained for the two types of disorder 
are qualitatively similar, the extension of regimes with 
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different spectral properties scales reciprocally with the 
disorder strength. 

While these results apply to the case of generic dis- 
order, there is one particular point in the phase dia- 
gram where drastically different behaviour is observed: 
In the combined limit half-filling/unitary scattering, the 
zero energy DoS exhibits a power law divergence H. In 
YAGHK the peculiar properties of this limit, henceforth 
referred to as the 'unitary limit' for brevity, were related 
to the formation of a second class of soft modes, existing 
in parallel to the particle-hole modes discussed above. 
It was shown that that perturbative inclusion of these 
modes leads to positive logarithmic corrections to the 
DoS. However, as with the generic particle- hole modes 
discussed above, no safe conclusions on the profile of the 
DoS can be drawn on the basis of pure perturbation the- 
ory. 

Below we will show that the emergence of a new class of 
soft modes is indicative of critical behaviour. The system 
at unitarity and e = sits on a critical line (parameter- 
ized through the ratio t/ A.) One has to emphasize that 
this fixed line lies outside the phase diagram of 'real' high 
T c materials; cuprates at half filling are anti-ferromagnets 
rather than superconductors. However, in the vicinity of 
the critical regime spectral and transport observables will 
exhibit anomalous scaling behaviour traces of which may 
(or may not) be observable in real systems. 

But what, then, is the origin of the exotic behaviour? 
In YAGHK the formation of the anomalous soft modes 
was traced back to the nesting symmetry of the half filled 
band. The 2x2 matrix Gorkov Hamiltonian function of 
the clean system obeys the symmetry 

AA:ff (k) = - J ff (k + Q), (9) 

where Q = (jr, n) is the nesting vector and a square lat- 
tice with unit lattice spacing is understood. Notice that 
what matters in the description of spectral/transport 
properties is not the Hamiltonian but the Green func- 
tion G — (e + fJ,o~3 — H)^ 1 (a~i, i = 1, 2, 3: Pauli matrices 
in particle-hole space), i.e. for finite energy arguments, e, 
or chemical potentials, the symmetry Af is effectively 
broken. 

Now, we are dealing with a system where translational 
invariance is broken by disorder, i.e. it seems more nat- 
ural to express the symmetry Af in real space. Noticing 
that translation by Q amounts to multiplication of every 
other site of the bipartite lattice by —1, we get 

Af:[Ho,af}+ = 0, (10) 

where [ , ]+ is the anti-commutator and the action of 
af is to multiply all sites belonging to one of the two 
sublattices A and B nested into the bipartite host lattice 
A\JB by -1. 

The representation ( |l(i| ) will not only turn out to be 
technically convenient, it also provides the clue as to 
the special role played by unitary scatterers. Evidently, 
a generic disorder potential superimposed on the clean 



Hamiltonian Hq will not respect the symmetry Af. With 
unitary, i.e. 'infinitely strong' scatterers, however, the 
situation is different. In a sense, such impurities burn 
local holes into the system. More technically, the po- 
tential v created by unitary scatterers is strong enough 
to locally force the amplitude of wave functions in the 
physical Hilbert space of the system to zero. Now, mul- 
tiplication of zero by ±1 does not change anything which 
is why the symmetry Af is effectively left intact by uni- 
tary scatterers. 

In which way does the symmetry Af affect the long 
range behaviour of the system? Temporarily forgetting 
about superconductivity, let us begin by briefly recapitu- 
lating the role of A/"-symmetries in conventional metallic 
systems. Building on early work by Wegner and Opper- 
mann ]l7j , Gade [18] realized that sublattice symmetries 
drastically affect the properties of disordered systems in 
the middle of the tight binding band. Specifically, (i) the 
middle of the band supports a delocalizcd phase and (ii) 

the DoS diverges as a power law v{e) ~ \e\~ 1 e~ K *v ln (l e l) j 
with some non-universal parameter k, upon approach- 
ing the band center. These phenomena reflect the fact 
that the e = system with A/"-symmetry belongs to a 
non-Wigner-Dyson symmetry class, viz. class ^4111 in the 
classification of Ref. Jl!| . Breaking this symmetry, finite 
values of the parameters e or /i induce a crossover to the 
Wigner-Dyson unitary class (or class A), i.e. the symme- 
try class of generic electronic systems with broken time 
reversal invariance. 

To obtain these results Gade described the normal 
metal sublattice system in terms of a (boson replica) ef- 
fective field theory, which was then subjected to a renor- 
malization group analysis. For future reference we note 
that this theory has a supersymmetric extension defined 
through Z = J VTexp(-S[T]) and 

S[T} = 

= IT / d2r [~ D str ( dTdT ~^ ~ 2ie str ( T + T ~ 1 )] + 
+ c Jd 2 r [st^rar- 1 )] 2 , (11) 

where c is a non-universal disorder-dependent constant 
and the fields T take values in the supergroup GL(1|1) 
(i.e. the group of invertible two-dimensional supermatri- 
ces.) 

While Gade obtained the field theory ( |TT| ) (or rather 
its boson replica analog) for the specific model of elec- 
tron hopping on a bipartite lattice, it is now understood 
that there is a more general way of looking at the prob- 
lem (cf., e.g., the discussion in Ref. pb]| ). Field theories 
like ([ll]) generally describe the long range behaviour of 
Hamiltonians of symmetry class AHI. A Hamiltonian be- 
longs to this symmetry class if it can be represented in a 
block form, 
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where Z are arbitrary matrices. (Indeed, the symmetry 
( |io| ) states that the Hamiltonian assumes a block off- 
diagonal form when represented in a sublattice basis.) 

The generalized interpretation helps to understand 
what happens with the sublattice system when the super- 
conductive aspects of the problem are switched on. In- 
deed, it is straightforward to show (cf. Appendix that 
(i) the superconductor sublattice Hamiltonian can be 
transformed to a representation ( fl2| ) and (ii) that there 
are no other symmetries in the problem. This implies 
that the superconductor with A/"-symmetry is again de- 
scribed by the low-energy effective action (|ll V a claim to 
be verified by explicit construction in section [V . There is 
one difference to the previously discussed case, viz. that 
the fields T now take values in the larger group GL(2|2) 
(i.e. extra degrees of freedom are needed to accommo- 
date the superconductor symmetries (Q)-) However, as 
far as the long range properties of the system are con- 
cerned, an extension of the field space to GL(n|n) (here 
n = 2) is inconsequential: As shown by renormalization 
group methods in Ref. [frsf , 



t> for zero excitation energy, e 
ical. 



0, the system is crit- 



> Upon approaching zero energy, the DoS diverges as 

Ky /-ln(\e\) 

i/(e) n . (13) 



where k is some non-universal, disorder depen- 
dent scale that cannot be reliably determined from 
the RG analysis. For sufficiently small energies, 
e < exp(— k 2 ), the DoS scales in a near power-law 
manner, v(e) ~ |e| _1 . 

o The conductance does not renormalize, i.e. the sys- 
tem is a (thermal) metal. 

t> All three, finite values of e, the chemical poten- 
tial [i, and deviations off unitarity, v^ 1 > 0, rep- 
resent relevant perturbations that drive the system 
towards the generic superconductor symmetry class 
CI. 

While the first three of the characteristics listed above 
can directly be inferred from Gade's analysis, the last 
needs some clarification. As shown in section IV, devia- 



tions off the unitary limit lead to the appearance of an 
extra operator 



s s b — r 



J d 2 r str (TT S + (TT S ) _1 ) 



(14) 



to be added to (pT|). Here, T = T(/i, v x ) is a symmetry 
breaking parameter whose explicit dependence on \x and 



v is given in Eqs. (|43|) and ( |46| ) below. Further, X s 
denotes a kind of generalized matrix transposition. More 
important than the explicit definition of this operation 



(see Eq. (g6[) below) is that finite values of T enforce 
T s = T -1 which is the defining relation of the orthosym- 
plectic group OSp(2|2) C GL(2|2). In other words, the 
action S s b induces a crossover between the extended sym- 
metry at the unitary limit (AEII) and the generic symme- 
try class CI When embedded into an RG analysis, 
the operator S s b represents a relevant perturbation with 
scaling dimension 2 (by elementary power counting.) I.e. 
deviations off the unitary limit rapidly increase under the 
RG operation and drive the system towards the phase of 
the generic superconductor. 




t-AIII 



FIG. 1. Qualitative diagram of regimes with different be- 
haviour of the DoS v(t). Horizontal axis: energy, e; verti- 
cal axis: A/"-breaking by finite chemical potentials /i and/or 
non-infinite scattering strengths v . On the abscissa the sys- 
tem is in class AIU, everywhere else in CI. For large values of 
e the DoS is well described by SCTA. The energy below which 
quantum corrections begin to play a role, increases with the 
anisotropy parameter A/t. Depending on whether the prod- 
uct evo is smaller or larger than the symmetry breaking pa- 
rameter F(fi,v~ 1 ) (the crossover region indicated by the solid 
line), quantum corrections to the DoS are negative or posi- 
tive. At the origin of the diagram (F, e) = (0,0), the system 
is critical and the DoS diverges. 



We finally discuss how these findings compare to the 
results of previous analyses, notably PL and YAGHK. 
To leading order, the divergence of the DoS, ([l3]), agrees 
with the result of PL. There are, however, some discrep- 
ancies on the level of logarithmic corrections to powerlaw 
scaling (PL find v(e) ~ l/|e|(ln(|e|) 2 .) whose origin we 
do not understand. The positive logarithmic corrections 
to the SCTA-DoS found in YAGHK are perturbative pre- 
cursors of the full powerlaw divergence. As a new result, 
we find that transport in the band center is metallic. 
It is tempting to interpret this metallic behaviour as a 
consequence of resonant tunneling transport in the nar- 
row impurity band created by unitary scattering centers. 
However, this picture is purely speculative. Finally, we 
disagree with the assertion of PL, that the divergence 
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survives generalization to finite values of fi. Using RG 
language, we rather find that [i acts like a (dynamically 
increasing) mass scale terminating the RG flow for small 
enough e. More specifically, one expects (cf. Fig. [I]) 
that for large values of e (i) the DoS will be given by the 
SCTA result. Perturbative quantum corrections to the 
large energy asymptotics (ii) scale logarithmically with e 
and inversely proportional to the anisotropy parameter 
tj A 3> 1. In an intermediate energy regime, RG sum- 
mation over all quantum interference contributions leads 
to the increasing DoS profile ([l3|). Finally, (iii) below 
a crossover scale set by voe/T(fi, v^ 1 ) = 1, the DoS fol- 
lows the decreasing profile of the superconductor of class 
CI. Non-monotonous behaviour of this type has indeed 
been observed numerically in YAGHK. Whether or not 
traces of the intermediate regime (ii) , or even of the scal- 
ing law (|l^), might be observable experimentally cannot 
be decided by the present approach. (Mentioning exper- 
iment, one should also notice that quasi-particle interac- 
tions act as an additional ingredient spoiling sublattice 
symmetries. I.e. the Hartree-Fock potential seen by the 
quasi-particles is local in real space and, therefore, vio- 
lates @.) 

This concludes our qualitative discussion. Starting 
from a microscopic lattice Hamiltonian, the subject of 
the next two sections will be the construction the field 
theories (||) and (|ll|). Basically, this will amount to an 
embedding of SCTA-summation techniques into a field 
integral approach. (Readers not interested in technicali- 
ties may skip these sections.) 



III. CONSTRUCTION OF THE FIELD THEORY 

Consider the quasi-particle Hamiltonian of a two- 
dimensional lattice d-wave superconductor, 



ij 

II], = £yC3 + Aijai, 

H$ , = V?8 ij a 3 + V*v a + V* L <T 1 , 

where £ y - = t tj - fx, 



(15) 



(16) 



is a spinor comprising a spin-up hole and a spin-down 
particle, and o~i are Pauli matrices acting in particle-hole 
space. The clean part of the Hamiltonian, H°, contains 
the chemical potential, fj,, a nearest neighbour hopping 
matrix element, tij, and the d-wave order parameter, Ay. 
The Fourier components of these operators are given by 



£(/c) = t(cos k x + cos k y ) — 
A(k) = A(cos k x — cos k y ), 



/i. 



where t and A define the strength of hopping and order 
parameter, respectively. 

The stochastic part of the Hamiltonian, H dls , accounts 
for two different types of disorder. First, the poten- 
tial generated by N randomly placed impurities is repre- 
sented through 



N 



If =5>(r*-R;), 



(18) 



where {Rj} are the coordinates of the impurity centers. 
Configurational averaging over Vf , 



N 

n h I **> ( 



will amount an independent, or Poissonian average over 
impurity coordinates. 

Second, the random operators V sa , a = 1, 3 model the 
presence of residual statistical fluctuations superimposed 
on the hopping and the order parameter amplitudes, re- 
spectively. For convenience, we assume these amplitudes 
to be Gaussian distributed with zero mean and variance, 



gSwS jf S aa 



(19) 



The motivation behind introducing Gaussian bond dis- 
order in addition to the Poisson impurity system is not 
only to make the modeling more 'realistic'. Below we 
will employ the Gaussian disorder average as a vehicle to 
introduce the soft fields T. Averaging over the impurity 
system will then be done in an after-step |E2). (Notice 
that for the purposes of the construction, the variance g 
can be assumed to be infinitesimally weak.) 

To prepare the construction of a field theory of this 
system, we change to an 'off-diagonal' representation, 
where the role of the Pauli matrix 173 is taken over by 
o 1 . This is done by a unitary transformation H — > 
exp(— ij<Ti)H exp(i^CTi), after which 

H°=ia 2 + Aa u 

H dis = {V? + V S2 )a 2 + V sl o-i. 

The new representation changes the appearance of the 
basic symmetries (Q), viz. 

C :H = -a 2 H T a 2 , 

T:H = a l H T a l . (20) 
We next introduce a supersymmetric partition function 



Z[J] = I D(V>,^)exp (iip 



e-H + J 



(21) 



from which disorder averaged matrix elements of the re- 
tarded Gorkov Green function can be obtained by differ- 
entiation with respect to elements of the source matrix 
J. In @, Ime > and 



(17) 



i/' b 
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is a four-component superfield whose complex com- 
muting (anticommuting) components V b ("0 f ) carry a 
particle- hole structure like in dlq). While convergence 
of the integral requires Tp h — 4> , the Grassmann fields 
?/r and t/j are independent. For transparency, the lattice 
index summation is suppressed in the notation. Simi- 
larly, the source field J will be set to zero in much of our 
further discussion. 

To fully account for both symmetries (|I]), it is neces- 
sary to subject the functional integral to a doubling, or 
'charge conjugation' operation. Using that 



the functional integral can be re-written as, 



Z[0} = 



D\I/exp ( 



-H 



where the eight-component fields 



on the non-directed links (ij) of the lattice and acting in 
bf and cc space. Finally, 



{Qa)i 



where N% is the set of four nearest neighbours of site i. 

We next subject the functional integral to a station- 
ary phase approximation. Temporarily neglecting both 
the energy increment e and the impurity potential one 
finds that the variational equations = are solved 
by the matrix-diagonal and spatially uniform configura- 
tion (QotQs) = (cr§ c , 0) . Here, n(g) is some function 
of the disorder strength which vanishes monotonously in 
the limit g \ 0. 

As usual with non-linear o-models, the diagonal con- 
figuration is but one representative of an entire manifold 
of solutions of the mean field equation. To identify this 
manifold, we imagine our diagonal configuration substi- 
tuted into the second line of ( p5| ) and explore what hap- 
pens when the 'I'-field is subjected to a transformation 
(matrix structure in ph-space), 



t 
a2ip T 



* = (tp, -ip T a 2 ® o- 



3 bf ) 



(22) 



Zi 



To 



To 



Here, the matrices <7i, of, and of c are Pauli matrices 
acting in particle-hole, boson-fermion, and in the newly 
introduced charge-conjugation space, respectively. The 
fields \& and ^ are related to each other by 



$ = $ a2®i", 



* = -02 ® 



where 



t = E\{ ® (ia?) + E\\ <g> of 



(23) 



(24) 



and is a matrix acting in bf-space and containing 
zeros everywhere save for a unity at position ij. Fortu- 
nately all we need to know about the ^'s is encapsulated 
by (|23J), i.e. there will be no need to explicitly unravel 
the entangled eight component structure of these fields. 

Continuing with the construction, we next average the 
functional integral over Gaussian disorder and decou- 
ple the resulting ^-interaction by means of a Hubbard- 
Stratonovich transformation: 



Z g [0] EE (Z[0]) f 



V(Q ,Q 3 )e 



--Y 

9 L^ij 



x / V^exp ( 



ea% c + 4Q + 4iQ 3 a 3 - £P 



= Jv(Q ,Q 3 )e-^ -W-^) x 
x exp -- str In (tof + 4Q + 4iQ 3 o 3 - ff p ) , (25) 

where H p = H L,_ . contains only Poisson disorder and 
Qa.ij, a = 0, 3 are four component supermatrices living 



Here, T\, . . . ,T 2 are spatially constant supermatrices of 
dimension four subject to the condition (to ensure com- 
patibility with (P3|) 



- x tT, 



Tf = Ti. 



(26) 



Remembering that H is off-diagonal in ph-space, we find 
that the transformation leaves the Hamiltonian invariant 
provided that 



Ti =T? = Tr 



(27) 



which is the defining relation of the supergroup OSp(2|2). 

We thus find that, notwithstanding the presence of the 
Poisson impurity potential in H p , the action has the con- 
tinuous group OSp(2|2)<g>OSp(2|2) as an invarian.ee man- 
ifold. While these transformations commute with the 
Hamiltonian, they change the diagonal saddle point con- 
figurations according to 



This equation identifies the set of all field configurations 
T ee r 2 ~VfTi e OSp(2|2) as the Goldstonemode man- 
ifold of the model. For fixed T 2 , the factor T 2 ~ of c S 
OSp(2|2) can be absorbed in Ti, i.e. the Goldstonemode 
manifold is isomorphic to a single copy of Osp(2|2) 
and can be parameterized by the representative T 
freely running through this group. Finally, disregard- 
ing massive fluctuations (i.e. field-configurations non- 
commutative with H) and using that the weight factor 
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str(Qg ii7 - + Q§ 



tion by T, we identify 



is invariant under transforma- 



j g = j T>T exp — — str In | < n , 



—H 21 i^T- 1 



as an effective functional capturing the low energy de- 
grees of freedom of the problem. 

We next have to face up to the analysis of the pointlikc 
disorder included in H p . More specifically, the object we 
need to compute is Z gp = (Z g ) p , expanded in slow fluc- 
tuations of the soft field T. Trading the average over a 
fixed number of impurities for a grand canonical average, 
it is possible to formulate an exact representation of the 
averaged functional |11|. However, transcendental in the 
fields T, the action of the grand canonical functional does 
not appear to be a particularly inviting starting point for 
subsequent low-energy expansions. 

Below we will adopt a different and less rigorous aver- 
aging scheme. The basic idea will be to employ elements 
of the standard approximations implied in the analysis of 
dilute systems of scattcrcrs, notably the self-consistent 
T-matrix approximation, in the expansion of the effec- 
tive action. To be more specific, what we want to do 
is expand the logarithm appearing in the action of Z g in 
(a) the 'mass' brought about by finite energy arguments e 
and (b) the finite action due to spatial fluctuations of the 
fields. The expansion of the logarithm in these pertur- 
bations, for a fixed realization of the impurity potential, 
leads to expressions like 

Z = fvT^J dfrMTin) J d 2 r 2 2 (T{v 2 ) 

where Oi(T) symbolically stands for a local operator de- 
pending on the fields T and a continuum approximation 
is understood. (Here, the attribute 'local' means local- 
ity on the scale of the elastic mean free path set by the 
scattering potential.) To make further progress with this 
expression, we rely on one central assumption: For a typ- 
ical configuration of coordinates (ri,r 2 , ■ . .), the spacings 
|rj — t*j I will be much larger than the correlation range 
£ of the impurity potential. For such configurations, the 
average can be assumed uncorrelated, i.e. 

(Oi(T(ri)0 2 (T(r 2 )...; r 

= (0 1 (T(r 1 )) p (0 2 (T(r 2 )) p . 

In a small subset of the integration volume, one or sev- 
eral integration variables come close to each other and 
averaging over impurity coordinates leads to residual cor- 
relations. This leads to the generation of composite op- 
erators, i.e. local operators of higher order in gradients 
and/or energies e. Relying on the fact that operators 
of this type will carry negative scaling dimension, we ne- 
glect statistical correlations and approximate the average 
by 

Zw f-DT ( d 2 r 1 (0 1 (T(r 1 )) p / d 2 r 2 (0 2 (T(v 2 )) p . . . , 



which amounts to replacing the average of the functional 
by an average of the action: 



VT x 



(28) 



x exp — - ( str In ( etTg 



inT 
-H p 

"21 



-H p 

■"12 

inT- 1 



We will re-assess the validity of this approximation be- 
low. 

The rest of the derivation is straightforward. We begin 
by computing the action S e corresponding to finite values 
of the energy parameter e. Expanding the action to first 
order in e and using that for sufficiently slow fluctuations 
the fields T and the Hamiltonian H p are approximately 
commutative, one obtains 



S £ = --str (TGu + T- x G 2 i) 



where G = (G%, 
G = 



in ~H P 2 
-H21 in 



(29) 



is the disordered Gorkov Green function and indices are 
in ph-space. 



» JL- 



+ 



b) 



FIG. 2. a) Diagrammatic representation of the SCTA series 
for the single particle Green function, b) The SCTA diffusion 
mode. 



Borrowing from previous analyses |g| , we compute the 
averaged superconductor Green function G by SCTA 
summation techniques ]2^] . The basic presumption un- 
derlying the SCTA (For a more comprehensive discus- 
sion of the approach and its application to the e?-wave 
superconductor see, e.g., Ref. p4|.) is that for decreas- 
ing impurity concentration, rii, correlations between scat- 
tering processes off different impurities do not affect the 
structure of the averaged Green function. Neglecting any 
amount of residual inter-impurity correlation, the SCTA 
approximates the self energy, S, of the Green function by 
the set of diagrams depicted in Fig. ^. Summation over 
these contributions leads to 

E = n,-T, 
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where the T-matrix (not to be confused with the slow 
fields T of the effective action approach) is defined 
through the Dyson equation 



T 



"IK7 2 + vga 2 T ■ 
f = (trVa 



Here, 5 = £ k G(k), 



G(k) 



£(k)a 2 - A(k)o-i - S 



(30) 



(31) 



and we anticipate that for a short range scattering poten- 
tial v(k) = v = const., the self energy is momentum inde- 
pendent. Expanding in Pauli matrices, £ = T,q + £202, 
and defining renormalized parameters k = k + i£ « i£ 
and £ = £ + £2, the integrated Green function g = 
.9o + 92 &2 assumes the form 



where 



.go = SqS, 92 = (m - s 2)3, 



(32) 



-5? 



£(k) 2 + A(k) s 



Linearizing the dispersion -ffo(k) around the four nodal 
points of the Fermi surface of the <i-wave superconductor 
one finds 



TrfA 



■In 



A 



(33) 



where A ^> £0 is a cutoff measuring the extent of the 
'linearization volume' in fc-space. With (pQ), and 



T = T + T 2 a 2 , 
we finally arrive at the self consistency equation 

go 



To 
T 2 



-So + - 92? 

v' 1 - 92 
-So + ~ 92? 



(34) 



where go, 2 depend on T$ through (|32|). There are two 
limiting cases, in which the solution of these equations is 
straightforward: In the Born limit, v^ 1 — g 2 3> |<?o|) 



„ Born 2 

To — > v g , 



m Born 

T 2 — > v, 



i.e. To is given by the familiar first order Born scatter- 
ing rate while T 2 is the averaged potential background. 
(Here we have neglected the chemical potential and the 
real part of the self energy for simplicity.) In the oppo- 
site, unitary limit, v^ 1 ~ g 2 — > 0, 



T, 



unitary 



T 2 un ^ y 0. 



With these results in place, we turn back to the analy- 
sis of S e . Expressing the trace over the average Green 
function through (B2J) and (J33|) one obtains 



(35) 



J d 2 r J d 2 k sir (T(r)Gu(k) + T- 1 (r)G 22 (k)) 
d 2 rstr (T + T- 1 ) , 



.■kevo 



where 



■ Im str (g) 



8|S | 
ttHA 



In 



A 

\Xo~\ 



(36) 



is the SCTA density of states at zero energy. 

We next consider the action Sr due to spatial fluctua- 
tions of the fields T. Turning back to (^) and temporar- 
ily setting t = 0we use the cyclic invariance of the trace, 
to re-write the action as 



S,=o[T] = \ 



str In 



-H\ 2 -T-\El 2 ,T\ 



"21 



(Notice that the impurity potential, diagonal in real 
space, commutes with the fields.) Due to [H® 2 ,T] — 
-iv ia (k) • VT + 0(d 2 T), where v 12 (k) = V k i?i 2 , the 
commutator between Hq and T contains derivatives of 
the slow fields and, therefore, counts as a small object. 
Expansion of the action to second order (The first order 
contribution vanishes by symmetry.) in the commutator 
leads to 



1 



str 



GaiCT- 1 via • VT) G 2 i(T _1 v 12 • VT) 



In the limit of negligibly small real space correlation of 
the potential (assumed throughout), v(k) = vo — const., 
vertex corrections to Sr vanish, i.e. 

S a = istr [G^T" 1 v 12 • VT) G 2 i(T _1 v 12 • VT)] . 

At this stage, the further evaluation of Sq seems straight- 
forward: One should substitute (|3l|), adopt the linear 
('nodal') approximation previously used in (|33"|), and do 
the momentum integral. However, there is a subtle prob- 
lem with this strategy. Straightforward substitution of 
the linearized dispersion leads to expression of the type 
• 00, where '0' is due to the /c-space symmetry of the 
integral and '00' due to ultraviolett divergencies. 

There are different ways to get around this problem. 
Abandoning the nodal approximation one may attempt 
to do the full momentum integral over the Brioullin zone. 
This integral is UV benign (We are working on a lattice!) 
but cumbersome. Alternatively, one may stick to the 
nodal approximation and implement one of several UV 
regularization schemes available for fermions with linear 
dispersion. Choosing the second route, and following a 
regularization procedure introduced in Ref. ||] , we obtain 







Sn[T] = - 



d 2 r str (OTdT- 1 ), 



(37) 



where the - disorder independent - coupling constant is 
related to the parameters t and A through (Q). Combi- 
nation of Eqs. ( p5[)_ and ( p7[ ) finally yields the effective 
low energy action(0) discussed in the introduction. 

At this stage, it is instructive to make explicit con- 
tact with the perturbative analysis of YAGHK. To this 
end, we expand the fields as T — exp(W) = 1 + W + 
W 2 /2 + . . ., where W € osp(2|2) are group generators of 
OSp(2|2). Expansion of the action to second order in W 
then obtains 



S[W] 



1 

2^ 



d 2 rstr (W (-Dd 2 - 2ie) W) + 



(38) 



where the kernel II -1 = —Dd 2 + 2ie is the inverse of 
the SCTA particle-hole diffusion mode depicted in Fig. 
H b) and derived diagrammatically in YAGHK. (That 
the quadratic action correctly describes diffusion on the 
background of the isolated impurity centers is evidence 
for the validity of the statistical assumptions underly- 
ing the derivation of the action.) To compute the first 
order quantum interference corrections to the DoS, one 
expands the functional expectation value (0) to second 
order in W and integrates. This produces the logarith- 
mically divergent results of YAGHK. For a list of more 
sophisticated schemes of evaluating the field theory, we 
refer back to the discussion of O. 



IV. THE HALF-FILLED BAND 

Building on the analysis of the previous section, we 
next consider the consider the unitary limit (by our abuse 
of language, the combined limit half filling/infinitely 
strong scattering.) To understand what will change, let 
us g o back to the level of the prototypical field theory 
( |25| ) and consider the bi-linear form 

Mo*- 

As discussed in section ||, the clean Hamiltonian Hq pos- 
sesses the symmetries C,T , and the sublattice symmetry 
Af. As we are aiming to describe a situation where not 
only C and T, but also Af is left intact, let us explore 
how this symmetry will change the internal structure of 
the field theory. 

As discussed in section [HI the symmetry Af is linked 
to the existence of two nested sublattices A and B. A 
superconductor Hamiltonian describing hopping between 
these lattices can be symbolically represented as 



Ho = hi ® o\ 



er 2 , 



(39) 



where the lattice Hamiltonians hi obey [^,03]+ = 0, 
% = 1,2, and af is defined in (JToj) . Likewise, the fields 
\& and ^ carry a two-component structure in sub-lattice 



space. The important point now is that for a theory glob- 
ally compatible with Af, the internal symmetry is larger 
than that discussed in the previous section. Indeed, it 
it straightforward to verify that ^S/Hq^S remains invariant 
under transformations 



8—1 



* = T#, (40) 



Tn 



TT 1 



*T. 



(41) 



where the matrix structure in sublattice and particle- 
hole space is defined through (Ap, Ah, Bp, Bh) and T 1>2 
are two independent matrices drawn from GL(2|2) D 
OSp(2|2). In other words, the symmetry of the theory 
has grown from OSp(2|2) x OSp(2|2) to the larger man- 
ifold GL(2|2) x GL(2|2). 

Having understood that, we proceed to explore (i) the 
structure of the fields entering the functional integral, 
(ii) the form of the low energy action, (iii) the role of 
finite chemical potentials /i, and, most importantly, (iv) 
why the symmetry Af survives generalization to unitary 
disorder. 

Temporarily ignoring the Poisson disorder, we notice 
that, as before, the diagonal configuration ina^ c solves 
the saddle point equation. Transformation by a smoothly 
fluctuating field of the structure (f40|), then leads to 



rps 



(42) 



which identifies T = Tx<jfT% £ GL(2|2) as the soft field 
of the theory and answers (i). 

Turning to (ii), we note that with the results of the 
previous section in place, the low energy structure of the 
action can be fixed by symmetry considerations, plus a 
few algebraic manipulations. Rather than explicitly ex- 
panding in slow fluctuations of the field T, we use that 
the low-energy action S[T] we are looking for must col- 
lapse to (^) if the symmetry Af is broken. This con- 
dition alone implies that the action of the system is 
given by (|TT|). Indeed, for fields drawn from the subset 
OSp(2|2) C GL(2|2) (i.e. the transformation group of the 
non-A/"-symmetric model) S[T] reduces to the action de- 
rived in the previous section. Further, S[T] contains all 
two-derivative operators compatible with the gross sym- 
metries of the system. The only coupling constant that 
cannot be fixed by simple symmetry reasoning is c. How- 
ever, for the purposes of the present analysis, the bare 
value of this constant is inessential. 

We next ask, (iii) , how finite values of the chemical po- 
tential 11 affect the picture. For finite /x, and still ignoring 
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the presence of the Poissonian disorder, the intermediate 
action assumes the form (cf. the fourth line in Eq. (I2F 



str In 



fieri - Hq 



One verifies that the isolated Pauli matrix o\ is incom- 
patible with the full set of transformations encapsulated 
in T; using the language of (|42]), only transformations 
with T s = T- 1 OSp(2|2) commute with o x . To 

quantitatively describe the symmetry crossover, one ex- 
pands the action to second order in the symmetry break- 
ing parameter \x. (The first order term vanishes after 
integration over the eigenvalues of Ho-) This obtains the 
symmetry breaking operator ([l4|), with coupling 



r = r a 



tA' 



(43) 



We finally, (iv), need to address the question of disor- 
der scattering. The generic Poisson potential V p ( |l8| ) is 
clearly incompatible with the lattice symmetry M . How- 
ever, for a singular potential, strong enough to locally 
eradicate all wave-function amplitudes in the low-energy 
sector of the Hilbert space, there is a chance for J\f to re- 
main effectively intact. An infinitely strong local poten- 
tial punches holes (sites where all wavefunctions vanish) 
into the system. Multiplication of these sites by a phase 
factor ±1 (the action of af) is a redundant operation, 
i.e. Af still holds. 

To quantitatively verify that this exceptional scenario 
is realized for unitary scattering at half filling, we con- 
sider the Green function operator 



We next employ the T-matrix equations ( p4[ ) and fia ) 
to explore the invariance properties of the theory un- 
der transformations T. Due to THoT^ 1 = Ho, the non- 
disordered Green function Go PI transforms as Go[T] = 
T -1 G f o[l]T -1 . Had the disordered Green function the 
same property, the theory would be invariant under 
the full set of transformations T. From Eq. (^5|) we 
see that invariance of the Green function is ensured if 
Ej[T] = T£[1]T. But if the Green function was invari- 
ant, we would have, 

E. t [T] = [fl-Va -T^GplJiiT- 1 ]" 1 
= T[ur 1 T<r 2 T-G[l]«]" 1 T. 

The last line shows that for generic values of the impu- 
rity strength, the self energy operator does not have the 
required invariance properties. Only in the limit of an 
infinitely strong potential Vi — > 00, the offending term 
T02T drops out and the action is invariant under the 
full set of transformations T. This is the unitary limit 
where the system is described by the soft action (p"l|). 
(Notice that for transformations drawn from the sub- 
group OSp(2|2) (i.e. the symmetry group of the non-A/"- 
symmetric superconductor) v~ 1 Ttj2 r F — 02, i-e. in- 
variance under the reduced set of transformations holds 
regardless of the impurity strength.) 

For scattering close to the unitary limit, v S> 1, an 
expansion of the action to second order in v" 1 obtains, 
again, the operator (p3), this time with coupling 



G\T] = inTa^T- 



H p 



as the central building block entering the expansion of 
the action. Notice that G[T] differs from the 'bare' Green 
function G, Eq. (^9|), in that the soft fields are included. 
The effect of impurity scattering on the local spectral 
weight is described by the self energy operator S[T] of 
the averaged Green function 



G[T} = 



ikTct§ c T - 



Ho -SfTl 



(44) 



To understand the connection between the impurity 
strength and the invariance properties of the Green func- 
tion, let us consider the self energy £j [T] due to a single 
local impurity potential Vi located at site i. (Effectively, 
this amounts to neglecting the effect of inter-impurity 
correlations on the wave-function amplitudes at individ- 
ual scattering sites.) Solution of the, now exact, Dyson 
series 



(v v) 2 



(46) 



The operator Jl4|), with coupling r = T\ + T 2 , describes 
the breaking ofAIII symmetry by finite chemical poten- 
tials and/or deviations off unitarity. An expansion of (|14| ) 
to second order in the field generators W (cf. Eq. (pq)) 
shows that for finite (/i, i; -1 ) the soft mode (the 'diffu- 
son') induced by the sublattice symmetry contains two 
mass terms, derived and discussed before by YAGHK. 
Whether or not traces of the AHI-type soft modes remain 
visible in the long range behaviour of transport and spec- 
tral observables depends on the relative value T/(i/oe). If 
this parameter is small, the mass of the superconductor 
soft modes is mainly set by deviations off zero energy and 
the (relatively small) additional breaking of the sublat- 
tice symmetry is inessential. In the opposite case, the 
sublattice symmetry is strongly broken and the system 
effectively in class CI. 



£j[T] = v, L a 2 + Via 2 G[T} i iE i [T} 



V. SUMMARY 



leads to 



£i[T] - [vT l a 2 -G[T\i$ 



This concludes our field-theoretical analysis of d-wave 
(45) superconductors with pointlike scatterers. We have 
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shown that in all but one point of a phase plane defined 
through impurity strength and filling factor, respectively, 
the system behaves similar to d-wave superconductors 
with generic, non-pointlike disorder. The DoS vanishes 
linearly upon approaching zero energy and eigenstates 
are localized (i.e. the system is a thermal insulator.) 

There exists one point in the phase diagram - infinitely 
strong disorder and half filling - where phenomenology of 
strikingly different type is realized: The DoS diverges as 
v(e) ~ |e| _1 /(ln |e|), where /(ln|e|) stands for logarith- 
mic corrections to power law scaling [Q , and zero-energy 
eigenstates are localized. Apart from deviations in the 
function / - whose origin remains obscure - the spectral 
profile found in the present analysis agrees with that of 
Ref. ||]. The identification of a metallic phase is a new 
result. 

We do not agree with the statement Q] that the pecu- 
liar behaviour realized at the critical point should be by 
and large insensitive to detuning of the parameters chem- 
ical potential/impurity strength. The question whether 
traces of critical behaviour - e.g. regions where the DoS 
increases upon lowering the energy - might be observable 
under 'realistic conditions' is not for the present approach 
to decide. However, given that the very formation of a 
superconducting phase requires a finite concentration of 
dopands (i.e. deviations off half-filling) , and that residual 
quasiparticle interactions spoil the sublattice symmetry, 
it seems likely that the 'real' superconductor will display 
behaviour characteristic for the generic symmetry class. 

Note added: Shortly before completion of this 
manuscript Ref. |^5| appeared. The authors of that refer- 
ence discuss structure and renormalization of supercon- 
ductor cr-models with sublattice symmetry (including the 
case of broken time reversal invariance.) Where the same 
symmetry classes are considered, the results of Ref. f25|| 
and the present paper agree. 

It is a pleasure to acknowledge discussions with C. 
Mudry, B. D. Simons, and M. Zirnbauer. 



APPENDIX A: BLOCK REPRESENTATION OF 
THE SUPERCONDUCTOR WITH NESTING 
SYMMETRY 

We show that the Hamiltonian H of the superconduc- 
tor with nesting symmetry can be transformed to a block 
off-diagonal representation. Due to the symmetry (|lC|), 
H can be represented through ([l2|), where the matrix 
structure is in sublattice space. The superconductor sym- 
metries (0) further imply that the block-matrices Z obey 



Z = Z* 



i.e. H has the form 



H 



OyZOy 



hi 
hi 



hi 
-h 2 



hi 
hi 



where the sub-structure of the blocks is in ph-space and 
the operators hi^ are real. A unitary transformation in 
ph-space brings H into the form 



H = 



Z* 



z 1 



where we have defined the complex operator Z = hi—ih^. 
Finally, changing the order of blocks, H can be rewritten 
as 



H 



z* 



z 1 



(Apart from an inessential transposition in the lower left 
hand block,) these are two identical copies of a class Alll- 
Hamiltonian. The replicated appearance of two AIII- 
Hamiltonians implies that our field theory takes values in 
the doubled field manifold GL(2|2) (instead of GL(1|1) 
for just one Hamiltonian.) 
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